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ABSTRACT 

As first Paper of a series devoted to study the old stellar population in clusters and fields in the Small 
Magellanic Cloud, we present deep observations of NGC 121 in the F555W and F814W filters, obtained 
with the Advanced Camera for Surveys on the Hubble Space Telescope. The resulting color-magnitude 
diagram reaches ~ 3.5 mag below the main-sequence turn-off; deeper than any previous data. We 
derive the age of NGC 121 using both absolute and relative age-dating methods. Fitting isochroncs 
in the ACS photometric system to the observed ridge line of NGC 121, gives ages of 11.8 ± 0.5 Gyr 
(Teramo), 11.2 ± 0.5 Gyr (Padova) and 10.5 ± 0.5 Gyr (Dartmouth). The cluster ridge line is best 
approximated by the a-enhanced Dartmouth isochroncs. Placing our relative ages on an absolute 
age scale, we find ages of 10.9 ± 0.5 Gyr (from the magnitude difference between the main-sequence 
turn-off and the horizontal branch) and 11.5 ±0.5 Gyr (from the absolute magnitude of the horizontal 
branch), respectively. These five different age determinations are all lower by 2 - 3 Gyr than the ages 
of the oldest Galactic globular clusters of comparable metallicity. Therefore we confirm the earlier 
finding that the oldest globular cluster in the Small Magellanic Cloud, NGC 121, is a few Gyr younger 
than its oldest counterparts in the Milky Way and in other nearby dwarf galaxies such as the Large 
Magellanic Cloud, Fornax, and Sagittarius. If it were accreted into the Galactic halo, NGC 121 would 
resemble the "young halo globulars", although it is not as young as the youngest globular clusters 
associated with the Sagittarius dwarf. The young age of NGC 121 could result from delayed cluster 
formation in the Small Magellanic Cloud or result from the random survival of only one example of 
an initially small number star clusters. 

Subject headings: star clusters: ages — star clusters: individual (NGC 121) — galaxies: Magellanic 
Clouds — galaxies: stellar content — stars: horizontal branch — stars: late-type 
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1. INTRODUCTION 

Characterizing old stellar populations provide impor- 
tant constraints on the early star formation histories of 
galaxies. Only the satellite galaxies of the Milky Way 
(MW) are sufficiently close to resolve individual stars 
well below the oldest main-sequence turn-offs, which 
is a pre-condition for accurate photometric age dating 
of old stellar populations. All Local Group galaxies, 
for which adequate data exist, appear to conta in stars 
older than 10 Gyr (jGrebel fc Gallagherl l20f)l . This 
result is based on main-sequence turn-off photometry 
of globular clusters and field populations in Galactic 
satellit es and a few more distant Local Gr oup galaxies 
(e.g., iBrown et all 120071 : ICole et al.|[2007D . as well as 
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the detection of horizontal branch stars (including RR 
Lyrae variables) i n the Local Group and beyond (e.g., 
i Held et all 120011: [Harbeck et alj 120011: ISaraiedini et all 
2002: IClementini et al.ll2003t IPritzl et al.H2004lh 

Globular clusters are preferred as the basis for old stel- 
lar population age tracers since they are usually single- 
age, single-metallicity objects facilitating comparative 
studies. Moreover, while g lobular cluster system s ex- 
hibit a range of ages (e.g., iDe Angeli et~afll2005r ). the 
oldest ones may belong to the most ancient surviving 
stellar systems to have comp leted their formati on in 
the youthful Universe (e.g., iMoore et alj [2006). In 
those nearby galaxies where relative age dating based 
on main-sequence photometry was carried out in com- 
parison to the oldest globular clusters in the Milky 
Way, no age diffe rence within the measurement accuracy 
was found (e.g.. jGrebel k, Galla gher 2004; B rown et al.1 
120071 ICole et al.ll2007l and references therein). The rel- 
ative age dating of the oldest identifiable Population II 
objects thus indicates a common epoch of substantial 
early star formation in the Milky Way and its compan- 
ions, although information about a putative, even older 
Population III remains to be uncovered in these objects. 

A galaxy that may not share this common epoch of 
early star formation - at least not with resp ect to its 
globular clusters (e.g., ISaraiedini et al.lll998h - is the 
Small Magellanic Cloud (SMC) 8 . The SMC is one of the 

8 There may be additional exceptions in more distant dwarf ir- 
regular galaxies regarding the common epoch of earliest Popula- 
tion II star formation, although also these galaxies evidently con- 
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closest and therefore best studied dwarf galaxies orbiting 
our Galaxy. 

While the SMC hosts a large number of intermediate- 
age and young star clusters, it only contains one "old" 
globular cluster, NGC 121, which is also the most massive 
star cluster. NGC 121 is located ~ 2.4°(~3 kpc) west 
of the SMC bar at (a J200 o.o, ^2000.0) = (0 ft 26 m 47.0 s , 
-71°32'12.0"). 

NGC 121 is the only cluster in the SMC that is suffi- 
ciently o ld to have developed an extended red horizontal 
branch (|Strvker et alj I1985T ) and to contain RR Lyrae 
stars. Indeed, whether or not to call a star cluster a 
globula r cluster is a matter of d efinition. In this case we 
refer to ISalaris fc GirardU (|2002f ) who consider Lindsay 1 
as having a stumpy red clump and not a red horizon- 
tal branch. Three RR Lyrae stars wer e disco vered in 
NGC 121 bv lThackeravl dl958h . iGrahaml (fl975h found a 
fourth RR Lyrae variable in the cluster and an additional 
75 in a 1 x 1.3 square degree field centered on NGC 121. 
Studies of various clusters in the LMC and in the MW 
showed that the presence of RR Lyrae variables indicates 
that the parent population is as old as or older than ~ 
10 Gyr. 

An important question is whether NGC 121 is as old 
as the typical old globular clusters in the Large Mag- 
ellanic Cloud (LMC) and in the MW. Previous stud- 
ies found ages rang i ng from 8 to 14 Gyr for NGC 121 
dStrvker et all Il98a iWalkeri H99ll iMiehell et all Il998l : 
lUdalskilll998l : iShara et al.lll998fc IDolphin et alJl200l us- 
ing a variety of different techniques. Studies based on the 
deepest available color-magnitude diagrams from Hubble 
Space Telescope (HST) observations with the Wide Field 
and Planetary Camera 2 (WFPC2) indicate an age of 10 
to 10.6 Gyr for NGC 121, suggesting that this globular 
cluster is several Gyr younger than the old est globulars 
in other nearby galaxies and in the MW (|Shara et al.l 
[TQM IDolphin et alJ[200lh . 

The capabilities of the Advanced Camera for Sur- 
veys (ACS) provide an improvement in both sensitiv- 
ity (depth) as well as angular resolution, which is essen- 
tial for a reliable photometric age determination in this 
dense star cluster. Here we present deep photometry of 
NGC 121 obtained with ACS aboard the HST. We de- 
termine the age of NG C 121 utilizing both ab solute and 
relative methods (e.g., IChabover et al.| [l996). The cur- 
rent study is the first in a series of papers based on HST 
studies of rich intermediate-age and old star clusters in 
the SMC. 

In addition to NGC 121, six intermediate-age SMC 
star clusters have been observed as part of our pro- 
gram: Lindsay 1, Kron3, NGC 339, NGC 416, Lindsay 38 
and NGC 419. We will derive fiducial ridgelines and fit 
isochrones to obtain accurate ages for each cluster using 
the same reduction techniques and isochrone models as 
described here (see §[2]|4|), and will present our results in 
future papers. In Table [1] we list the cluster identifica- 
tion, date of observation, passband, exposure times and 
location of all clusters in our HST program (GO-10396; 
principal investigator: J. S. Gallagher). 

In the next Section we describe the data reduction pro- 
cedure. In § [3] we present the color-magnitude diagram 
(CMD) of NGC 121 and discuss its main features. In §H 

tain old populations (e.g., IGrebeUl200ll : IMakarova e t al. 2002) 



we describe our age derivation methods and present our 
results. 

2. OBSERVATIONS AND REDUCTIONS 

The SMC cluster NGC 121 was observed with HST's 
ACS on 2006 March 21 as part of our program focused 
on star clusters and field stellar populations in the SMC. 
The program aims at exploring the star formation his- 
tory and properties of the SMC using both a number of 
carefully selected clusters and field regions. For NGC 121 
we obtained imaging in the F555W and F814W filters, 
which resemble the J ohnson V and I filte rs in their pho- 
tometric properties (jSirianni et al. 2005). The images 
were obtained using the Wide Field Channel (WFC) of 
ACS and cover an area of 200" x 200" with a pixel scale 
of ~ 0.05 arcsec. One set of exposures was taken at the 
nominal position of the cluster center. Four long expo- 
sures were obtained in each filter for hot pixel removal 
and to fill the gap between the two halves of the 4096 
x 4096 pixel detector. Each pointing has an exposure 
time of 496 s in the F555W, and 474 s in the F814W 
filter. Moreover, two short exposures were taken in each 
filter with an exposure time of 10 s in F555W and 20 s 
in F814W. 

The data set was processed adopting the stan- 
dard Space Telescope Science Insitute ACS calibration 
pipeline (CALACS) to subtract the bias level and to ap- 
ply the flat field correction. For each filter, the short 
and long exposures were co-ad ded independently using 
the MULTIDRIZZLE package (|Koekemoer et alJl2002[ ). 
With this package the cosmic rays and hot pixels were 
removed and a correction for geometrical distortion was 
applied. The resulting NGC 121 data consist of one 40 s 
and one 1940 s exposure in F555W and one 20 s as well as 
one 1896 s exposure in F814W. The two short exposures 
allowed us to measure brighter stars that are saturated 
in the long exposures. 

The photometric reductions were carried out using 
the DAOPHOT package in the IRAF environment 9 . 
We discarded saturated foreground stars and back- 
ground galaxies using the Source Extractor package 
(|Bertin fc Arnoutslll996l K 

Due to the different crowding and signal-to-noise ra- 
tio properties of the long and the short exposure images, 
photometry involving point spread function (PSF) fitting 
was only performed on the long exposures. For the short 
exposures we used aperture photometry, which turned 
out to yield smaller formal errors than PSF photometry. 
We ran DAOPHOT on our data and set the detection 
threshold at 1 a above the local background level in or- 
der to detect even the faintest sources. The list of stars 
detected in the F814W image was then used as coordi- 
nate input list to identify the stars in the F555W im- 
age and serve as our coordinate reference. 49493 sources 
were found to be common to both long exposure frames. 
For these sources, we performed aperture photometry us- 
ing an aperture radius of 3 pixels. We then constructed 
a PSF by combining 150 bright and isolated stars that 
were distributed fairly uniformlly across the image. Fi- 

9 IRAF is written and supported by the IRAF programming 
group at the National Optical Astronomy Observatories (NOAO) 
in Tuscon, Arizona. NOAO is operated by the Association of Uni- 
versities for Research in Astronomy, Inc. under cooperative agree- 
ment with the National Science Foundation. 
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TABLE 1 

Journal of Observation 



Image Name Date Total Exposure Time R.A. Dec. 

Cluster (yy/mm/dd) Filter (s) 



NGC 121 



Lindsay 1 



Kron 3 



NGC 339 



NGC 416 



Lindsay 38 



NGC 419 



J96106030 
J96106040 
J96106010 
J96106020 
J96105030 
J96105040 
J96105010 
J96105020 
J96107030 
J96107040 
J96107010 
J96107020 
J96104030 
J96104040 
J96104010 
J96104020 
J96121030 
J96121040 
J96121010 
J96121020 
J96102030 
J96102040 
J96102010 
J96102020 
J96103030 
J96103040 
J96103010 
J96103020 



2006/03/21 



2005/08/21 



2006/01/17 



2005/11/28 



2006/03/08 



2005/08/18 



2006/01/05 



F555W 
F814W 
F555W 
F814W 
F555W 
F814W 
F555W 
F814W 
F555W 
F814W 
F555W 
F814W 
F555W 
F814W 



40.0 
1984.0 

20.0 
1896.0 

40.0 
1984.0 

20.0 
1896.0 

40.0 
1984.0 

20.0 
1896.0 

40.0 
1984.0 

20.0 
1896.0 

40.0 
1984.0 

20.0 
1896.0 
40.0 

1940.0 
20.0 

1852.0 
40.0 

1984.0 
20.0 

1896.0 



h 26 m 42.98 s 
h 26 m 43.26 s 
h 26 m 42.98 s 
/l 26 m 43.26 s 
h 03 m 53.19 s 
ft 03 m 52.66 s 
h 03 m 53.19 s 
h 03 m 52.66 s 
h 24 m 41.64 s 
h 24 m 41.92 s 
h 24 m 41.64 s 
h 24 ra 41.92 s 
h 57 m 47.40 s 
h 57 m 47.13 s 
h 57 m 47.40 s 
h 57 m 47.13 s 
l h 07 m 53.59 s 
l h 07 m 54.09 s 
l h 07 m 53.59 s 
l h 07 m 54.09 s 
,l 48 m 57.14 s 
h 48 m 56.76 s 
,l 48 m 57.14 s 
ft 48 m 56.76 s 
l h 08 m 12.53 s 
l' l 08 m 12.71 s 
l' l 08 m 12.53 s 
l h 08 m 12.71 s 



-71°32'16.54" 

-71°32'14.61" 

-71°32'16.54" 

-71°32'14.61" 

-73°28'15.74" 

-73°28'16.47" 

-73°28'15.74" 

-73°28'16.47" 

-72°47'47.49" 

-72°47'45.49" 

-72°47'47.49" 

-72°47' 45.49" 

-74°28'26.25" 

-74° 28' 24. 16" 

-74°28'26.25" 

-74°28' 24.16" 

-72°21'02.47" 

-72°21'01.79" 

-72°21'02.47" 

-72°21'01.79" 

-69°52'01.766" 

-69°52'03.07" 

-69°52'01.76" 

-69°52'03.07" 

-72°53 / 17.72" 

-72°53 / 15.49" 

-72°53'17.72" 

-72°53'15.49" 



nally, PSF photometry was carried out. 

The photometric calibration was accomplished by con- 
verting the magnitudes of the individual stars to the 
standard ACS magnitude system by using an aperture 
with a radius of 0.5" (or 10 pixels on the image), in com- 
bination with the aperture correction from the 0.5" aper- 
ture radius to infi nity and the synth etic zero points for 
the ACS/WFC (|Sirianni et al.ll200l . The aperture cor- 
rection was derived for each frame independently. The 
objects found in both images were cross-identified and 
merged with a software package written at the Bologna 
Observatory by P. Montegriffo (private communication). 
Altogether we were able to cover a luminosity range of ~ 
10 magnitudes after combining the resultant photometry 
of the short and long exposures. 

In Figure [T] we show the photometric errors assigned 
by DAOPHOT. For stars measured on the short ex- 
posures, the formal photometric errors remain negligi- 
ble over a wide range of magnitudes. In the long ex- 
posures, all the brighter stars with 771555 < 19-3 mag 
and m§i4 < 19.5 mag are saturated. At 777555. 814 ~ 
19.6 mag, the short exposure (blue dots) and aperture 
photometry from the long exposure (black dots) samples 
were combined, and for stars fainter than 777555,814 ~ 
22.2 mag, long exposure PSF photometry (red dots) was 
used. We chose where to cut between the aperture and 
PSF photometry catalogues based on the 777555 data and 
adopted the same value for msu so as to avoid a color 
slope associated with this division. For our study, we 
rejected all stars with a a error larger than 0.2 mag and 
a DAOPHOT sharpness parameter -0.2 < s < 0.2 in 
both filters. To obtain a superior CMD, we discarded all 



stars within a radius of 35" from the cluster center, which 
excludes the very dense core of the cluster. With this se- 
lection, our final sample contains 17464 stars common in 
both filters. 

3. THE COLOR-MAGNITUDE DIAGRAM 

The resulting color-magnitude diagram (CMD) of 
NGC 121 and its surroundings is shown in Figured Our 
CMD for NGC 121 reaches -3.5 mag below the MSTO 
(~ 0.5 magnitudes deeper than the previous deepest 
available photometry), which allows us to carry out the 
most accurate age measurements obtained so far. The 
CMD shows a well-populated main sequence (MS), sub- 
giant branch (SGB), red giant branch (RGB), horizontal 
branch (HB), and asymptotic giant branch (AGB). The 
gap on the RGB at 777555 ~ 20 mag is an artificial feature 
due to small number statistics resulting from our exclu- 
sion of crowded stars in the cluster center. NGC 121 
appears to be a single-age population object just as one 
would expect for a canonical star cluster. As expected, 
there is no obvious evidence for field star contamination 
by younger populations due to the location of NGC 121 in 
a low-density area in the outer parts of the SMC. Within 
the field of view of the ACS and at the high Galactic 
latitude of the S MC, Galactic foreground cont amination 
is very low (e.g., iRatnatunga &; Ba hcall 1985]). 

Another possible contamination source is the massive 
and extended Galactic globular cluster 47 Tu c, which has 
a tidal radius of 42.86 arcmin ([Harris! fl996h and an an- 
gular distance from NGC 121 of ~ 32 arcmin. 

We visually estimated the location of the center of 
NGC 121 on the image and selected all stars within an 
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Fig. 1. — Photometric errors assigned by DAOPHOT to stars 
in the short exposures (blue dots), in the aperture photometry 
from the long (red dots), and on the PSF photometry from the 
long (black dots) exposures. Note the very small formal errors in 
the aperture photometry of the short exposures. Stars brighter 
than ~ 19.3 mag in the long F555W exposure and brighter than 
~ 19.5 mag in the long F814W exposure are saturated and are 
therefore not shown. The lower envelope of the error distribution 
of the stars in the short and long F555W exposure (aperture pho- 
tometry) cross over at m555,8l4 = 19.6 mag, and in the aperture 
and PSF photometry at r7i555^i4 = 22.2 mag (also indicated by 
two thin vertical lines). Here the photometry of the short and long 
exposures was combined. For the F814W exposures we chose the 
same magnitude value in order to avoid introducing a color slope 
in the color-magnitude diagram of the resultant data set. 




m 555 _m 814 

Fig. 2. — Color-magnitude diagram of NGC 121 and its surround- 
ings. Stars within a radius of 35" from the cluster center have been 
discarded. All stars with "good" photometry (cr < 0.2 mag and 
0.2 > sharpness > —0.2) are shown; 17464 stars in total. Repre- 
sentative errorbars (based on the errors assigned by DAOPHOT) 
are shown on the left for the msss — msi4 color. 

annulus of 35" and 45" to create a bona fide sample. 
This CMD is displayed in Figure [3] There is no evidence 
for a binary sequence in NGC 121, but we cannot ex- 
clude their presence, due to the photometric error. The 
aforementioned traces of minor field contamination have 
mainly vanished. Due to crowding, incompleteness be- 




FlG. 3. — Color-magnitude diagram of all stars within an an- 
nulus between 35" and 45" of NGC 121. We used this CMD for 
the determination of a representative color-magnitude ridgeline of 
NGC 121 (cyan line). This CMD contains 5112 stars. Only stars 
with good photometry (cr < 0.2 mag and 0.2 > sharpness > —0.2) 
are shown. 

comes significant at the faint end of the MS: This af- 
fects particularly faint stars in the cluster center. Hence 
in Figure [3] the MS becomes less densely populated at 
fainter magnitudes. 

The red HB is well populated and extends into the 
RR Lyrae instability strip (Clementini et al. 2007, in 
prep). The presence of a red HB provides a cir- 
cumstantial suggestion that NGC 121 may be younger 
than old Galactic and LMC globular clusters, since 
the HBs of the oldest globular clusters tend to ex- 
tend farther into the blue (e.g.. | Qlszewski et alJH996t 
IQlsen et alJH99cl IMackev fc Gilmord 120041) . Red HBs, 
however, can also be due to a "second parameter" 
other than age affecting the HB mo r phology (e.g., 
Lee et al.l 119941: iBuonanno et alj Il997t lHarbeck et all 
20011 : ICatelan et al.l l200lf T Since a true HB is present, 
an age measurement for NGC 121 can be made using 
the AV^P age measurement, which we will do in § 14.21 
This method requires the determination of the appar- 
ent mean magnitude of the HB. Our data yield 771555^3 
= 19.71 ±0.03 mag for this observable, which is in 
agreement with the magnitudes determined by Shara et 
al. (1998), Alves & Sarajedini (1999) and Dolphin et 
al. (2001). 

At m S55 = 19.58 ± 0.03 mag we find the NGC 121 
RGB bump (771555. Bumv) which is 0.06 mag brighte r than 
the mag nitude found by I Alves fc Sarajed ini (1999). The 
difference in luminosity is due to the exclusion of the 
center stars. If we determine the m^ t Bump on the en- 
tire sample, we obtain 771.555 = 19.52 ± 0.04 mag, which 
is in excellent agre e ment with the magnitude found by 
lAlves fc Saraiedml (|1999t ). This feature is predicted by 
stellar evolution models, which also show that the lumi- 
nosity of the RGB bump is dependent on the metallicity 
and age of the cluster. When the metallicity is known 
the difference between Vhb and Vsump can be used as 
an age indicator. 

Above the MS turn-off, IShara et al] (|1998[ ) found 42 
candidate blue straggler stars (BSS). Evolved descen- 
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Fig. 4.— The CMD of NGC 121 with the best-fitting isochrones 
of two different models: The blue solid line shows the best-fitting 
Padova (Girardi, "private communication", Girardi et al. 2000) 
isochrone that is closest to the spectroscopically measured metal- 
licity of the cluster. The r ed solid line is the best-fitting Ter- 
amo (Pietrinfcrni ct al. 2004) isochrone approximating the known 
metallicity. Neither model is a-enhanced. The cyan solid line is 
our fiducial ridgeline. The fitting parameters are listed in the plot 
legend. 

dants of the BSSs are important as pos sible sources o f 
stars lying above the traditional HB ( e.g., Catclan 20051) . 
In our ACS study we recovered the Sh ara et alj (|1998f ) 
BSS sample and also found more stars in the BSS region 
of which some (about 20) turned out to be pulsating vari- 
ables (dwarf Cepheids). These stars will be discussed in 
Clementini et al. (2007, in preparation) where we will 
present the results of an HST study of variable stars in 
NGC 121. 
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4. AGE OF NGC 121 
Age Based on Isochrone Fits 



Age determinations of star clusters using isochrones 
depend crucially on the interstellar extinction, distance, 
and metallicity of the cluster, as well as on the chosen 
isochrone models. In fitting isochrones to the CMD of 
NGC 121 we adopted the spectroscopic metallicity mea- 
surement of [Fe/H] = -1.46 ± 0.10 from Da Costa & 
Hatzidimitriou (1998, see also Joh nson et al. 2004) on 
the metallicity scale introduced by IZinn &: Westl (|1984l ) 
(ZW84). The distance and the extinction were treated as 
free parameters. The SMC dista nce modulus is (m — M)p 
= 18.88±0.1 mag (60 kpc) (e.g., iStorm et al.ll2004h . but 
due to the large depth extension of the SMC along the 
line of sight we adjusted the distance modulus (m — M) 
to produce the best isochrone fits to our CMD data. 

For easier comparison to the isochrones, we first de- 
rived a fiducial ridgeline (Table [5]), which reproduces the 
mean location of the stellar distribution in the CMD (ex- 
empting the HB) . In order to determine the ridgeline, we 
separated the cluster center CMD into three sections: the 
MS, the SGB and the RGB. On the MS we determined 
the mode of the color distribution in magnitude bins of 
0.3 mag width. For the SGB, we performed a linear least 
squares fit of a polynomial of 5th order to a Hess diagram 
of this region in the CMD. Finally, the RGB was fit by a 



Fig. 5.— The N GC 121 CMD with the best-fitting Dartmouth 
IIDotter et al.ll2007l) isochrones overplotted in red. As before, the 
cyan line represents our fiducial for NGC 121. The fit parameters 
are listed in the plot. Note the excellent agreement of this a- 
enhanced isochrone with the observed CMD. 

third-order polynomial of the mean color, again in mag- 
nitude bins with a size of 0.3 mag each. The resulting 
ridgeline is shown in Fig. [3] as a cyan line. 

We fitted our 771555 vs. 771555-771814 CMD with three 
different isochrone models: Padova isochrones (Girardi, 
"private communic ation", Girardi et al. 2000) 10 , 
Teramo isochrones (|Pietrinferni et aill2004D . both with 
scaled solar isochrones ([a/Fe ] = 0.0), and Dartmouth 
isochrones (|Dotter et al.l l2007f ) with both [a/Fe] = 0.0 
and +0.2. The Padova isochrone grid has an age reso- 
lution of log(t)=0.05, the Teramo isochrone grid of 0.1 
Myr and the Dartmouth isochrone grid of 0.5 Gyr. Our 
adopted spectroscopic metallicity of [Fe/H] = —1.46 cor- 
responds most closely to Z = 0.0004 in the Padova mod- 
els, to Z = 0.0006 in the Teramo models, and to [Fe/H] 
= —1.49 in the Dartmouth models. All three sets of 
isochrone models are available in the standard ACS color 
system. 

We fitted a large number of isochrones using differ- 
ent combinations of reddening, age, and distance. For 
each set of models, we selected the isochrone that best 
matched the observed data (Fig. [4] Fig.[5|). 

First we discuss Figure _U Our best-fit age us- 
ing Padova isochrones is 11.2 Gyr with (m — M)o = 
19.05 mag and Ey-i — 0.05. The best fitting Teramo 
isochrone yields an age of t = 11.8 Gyr, (m — M)o = 
19.02 mag, and £V_/ = 0.04. On the MS, both the 
Teramo isochrone and the Padova isochrones trace the 
ridge line almost perfectly. At the faint end of the MS, 
the Padova isochrone continues further to the blue than 
the Teramo isochrone and our derived ridge line; how- 
ever, this only becomes more apparent at magnitudes of 
m 555 ^25.5 mag and below. Both isochones also provide 
an excellent approximation to the SGB and to the lower 
RGB up to about half a magnitude below the HB. 

At brighter magnitudes, the two isochrones deviate in- 
creasingly to the blue of the observed upper RGB. Here 
the Padova isochrone shows the strongest difference, de- 

10 http://pleiadi.pd.astro.it/isoc_photsys.02/isoc_acs_wfc/index.html 
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viating by approximately 0.38 mag in color from the ob- 
served tip of the RGB. The isochrone shows a magnitude 
for the base of the red HB that is about 0.5 mag fainter 
than the observed one. Unlike Teramo and Dartmouth, 
the Padova isochrone also models the AGB and its tip, 
which is ~ 1 mag brighter than the tip of the RGB. The 
Teramo isochrone is too blue by about 0.23 mag at the 
magnitude of the tip of the RGB and indicates a mag- 
nitude for the base of the red HB that is 0.2 mag too 
bright. 

If we had no prior knowledge of the metallicity of 
NGC 121 and were to use the upper RGB as a metallic- 
ity indicator, a better fit would be obtained by choosing 
isochrones of a different metallicity or a abundance. The 
problems of various isochrone models of given metallici- 
ties in reproducing the upper red giant branches of globu- 
lar clusters with the same metallicit ies are a well-known 
problem (e.g., iGrebell Il997l . 1 19991) . Our Figure H re- 
flect the general failure of the chosen stellar evolutionary 
models to simultaneously rep roduce the major features 
of CMDs (|Gallart et al.ll2005f) in spite of the excellent fit 
to the lower RGB, SGB, and MS. Fortunately the latter 
are the most age-sensitive features of the CMP. 

The isochrone model provided by iDotter et al.1 (|2007f ) 
with [a/Fe] = +0.2 yield the best fit to the CMD (Fig.©. 
The best-fit isochrone has the parameters t = 10.5 Gyr, 
(m — M)o = 18.96 mag, and Ey-i = 0.035, using the 
a-enhanced isochrones of [a/Fc] = +0.2. All the major 
features of the CMD are very well reproduced, including 
the upper RGB where the isochrone is offset slightly to 
the blue relative to the fiducial ridgeline. This offset is 
no more than 0.01 to 0.02 on average along the entire 
upper RGB; i.e., even the slope of the RGB is very well 
reproduced along its entire extent. Unfortunately the 
stellar evolution models used here terminate at the He 
flash, and therefore do not fit the HB or the AGB. 

Is our use of a-enhanced models justified? For 
NGC 121 a value of [Ca/Fe] = +0.24 has been mea- 
sured, which is similar to the outer LMC cluster Hodge 
11 and to the old Galactic outer halo clusters with 
[Ca/Fe] = +0.3 (| Johnson et al.ll2004f) . Consequently, we 
assume that NGC 121 also is enhanced in a-elements. We 
note that in this respect NGC 121 differs from the gen- 
eral trend observed in red giant stars in the LMC and 
in dwarf spheroidal galaxies, where the [a/Fe] ratios at 



Finally, we note that the derived reddeni ngs agree with 
the ex tinction Ay = 0.1 ± 0.03 fro m the ISchlegel et al.l 
maps. The reddening law of lO'Donnelll (|1994[ ) is 



Shetrone et al 



Johnson et al.l 



. iHill et all 


2000|; 


iPritzl et all 


2005; 



a dex than in the Galactic halo (e.p 
[200lt lFulbrighdr2 002 
20061 : iKoch et al.ll2007l) . 
When we adopt the values for distance and redden- 
ing, but fit the cluster with an isochrone scaled so- 
lar, NGC 121 gets a slightly older age of 11 Gyr. The 
isochrone model with [a/Fe] = 0.0 still provides a better 
fit than the Teramo or Padova models, but has an offset 
of ~ 0.05 on average along the upper RGB. Past studies 
found that a-enhanced models imply a higher luminosity 
and temperature for the same mass than the solar scaled 
models and therefore an older age fo r the same magni- 
tude (e.g., IVandenBerg et al.ll2000T ). The Dartmouth 
models show exactly the opposite behavior. This is be- 
cause in these models an increase in [a/Fe] is accompa- 
nied by a corresponding increase of the total metallicity 
Z, which makes the isochrones cooler at constant age and 
[Fe/H] . 



assumed. 

In the Figures [6] and [7] we show a range of isochrones 
for the three sets of stellar evolution models in order 
to illustrate the age uncertainty in a given model. The 
finally chosen, "best" isochrone is always displayed 
along with two younger and two older isochrones. For 
the cases of the Teramo and Padova models, the two 
isochrones that are one age step younger or older than 
the chosen, central isochrone provide an upper or lower 
envelope for the location of the high-density part of the 
SGB, the MS turn-off, and the base of the RGB. For the 
Dartmouth isochrones the outermost isochrones provide 
this envelope. Considering the high quality of the fit of 
the central isochrone in this CMD region in all models 
and the larger deviations of the adjacent isochrones, 
we estimate that the resultant age uncertainty is of the 
order of approximately ±0.5 Gyr for the Teramo and 
Dartmouth isochrones and ±0.7 Gyr for the Padova 
isochrones. 

4.2. Empirical Age Estimates 
4.2.1. Vertical Method 

To check the reliability of the isochrone ages, we 
use a reddening-independent method to derive relative 
ages of NGC 121. This method is also independent of 
the photometric zeropoint of our data. This "vertical 
method" relies on the fact that the absolute magnitude 
of the MSTO depends o n the age of the cluster (e.g., 
lAlves fc Sarajedimlll999| ). while the absolute magnitude 
of the HB remains approximately ag e-independent for 
clust ers older than t > 10 Gyr (e.g., iGirardi fc Salarisl 
I2001T ). We measure the apparent magnitudes of the 
MSTO and HB, i.e., Vto an d Vhb, in order to obtain 
the magnitude difference AV^q . With increasing age, a 
cluster has generally larger values of this parameter since 
the MS moves to fainter magnitudes. Unfortunately, the 
determination of these two points comes with significant 
uncertainties. AVrpQ is hard to measure accurately both 
because of the width of the HB in luminosity and the 
MSTO's vertical extent in the turn-off region. 

First, we calculate the magnitude difference Amfo 555 
bet ween the HB and the M STO, as originally described 
by llben fe Faulknerl (|1968f ). Because W555 is propor- 
tional to V and we are only interested in the magni- 
tude difference, which we measure at constant color, 
there is no need to transform the magnit udes from the 
ACS system to V and I magnitudes (|Sirianni et al.l 
l2005h . We follow the general definition of the MSTO 
as the bluest point along the MS; in our case repre- 
sented by the bluest point on the ridgeline. We find the 
MSTO at rriTo.,555 = 22.98 ±0.05 mag, (to 555 -m S u) = 
0.59 ± 0.005 mag. As described in Section [31 the mean 
HB magnitude is m,HB,555 — 19.71 ± 0.03 mag. This 
yields a Am^ 555 = 3.27 ±0.06 mag. Our result is only 
slightly lower than former values, e.g., A V^P = 3.32 
determ ined by[Strykcr et al. (19 8a), 3 -33 by Shar a et all 
(|1998Q or 3.29 bv lDolphin et al.l (|2001h . iBuonanno et~aTl 
(1989) published a mean value of AV^P = 3.55 mag for 
old Galactic halo globular clusters (GC). The Ato^q 555 
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Fig. 6. — The color-magnitude diagram of NGC 121 after zoom- 
ing in on the region of the main-sequence turn-off, subgiant branch, 
and lower red giant branch. In the upper panel, we show Ter- 
amo isochrones as solid lines, covering an age range of 10, 10.9, 
11.8, 12.6, and 13.5 Gyr. These are the age steps in which these 
isochrones are provided ( Pictrinfcrni ct al. 2004). The central 
isochrone is our chosen best-fitting isochrone. In the lower panel 
we show the same plot for Padova isochrones (solid lines) in the 
Padova age steps of 8.9, 10, 11.2, 12.6, and 14 Gyr (Girardi, "pri- 
vate communication", Girardi et al. 2000). All other parameters 
are the same as in Figs. l4l and |51 

of NGC 121 is 0.28 mag smaller than this value, which 
indicates that NGC 121 is younger than most of the older 
Galactic GCs. 

Walker (1992) (see also Buonanno et al. 1989) found a 
relation between age and Al^f based on a study of 41 
Galactic globular clusters. Adopting again a metallicity 
of [Fe/H] = -1.46±0.10, we find with the formula logi = 
-0.045 [Fe/H] + 0.37AVgf -0.24 an age of 10.9±0.5 Gyr 
for NGC 121. 



4.2.2. Age Estimate using Mv{HB) 

In order to co mpare our data directly with 
iDolphin et all (120011). we u se the age calibration provided 
bv lChabover et al r i[l996bfl as was done by these authors. 
For this purpose we must adopt the metallic i ty Fe /H] 
= -1.19±0.12 of lDa Costa fc Hatzidimitrioul (11998D on 
the m etallicity scale introduced bv ICarretta &: Grattonl 
1997fl (CG97) to be consistent with IDolphin et~aT1 
20011 ) 's calculation. Note that elsewhere in the pa- 



per we are using metallicities on the ZW84 scale which 
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TABLE 3 

Ages for NGC 121 derived in this paper 



24 



Fig. 7. — Same as Fig. [6] but for the ct-enhanced Dartmouth 
isochrones covering an age range of 9.5, 10, 10.5, 11, and 11.5 Gyr 
flDotter et al.|[200W) . 

agrees reasonably well with the spectroscop ically derived 
chemical abundances ( Joh nson et al.ll2004f ). Using this 
method and the higher metallicity we determine an age 
of 9.7 ± 1.0 Gyr for NGC 121 , which is similar to the 
age that iDolphin et al.l (|2001l ) found. We have to em- 
phasize that this age is younger than the absolute ages 
determined in Section [4. II due to the different metallicity 
scale. If we use our preferred ZW84 metallicity scale with 
this method we obtain an older age of 10.8 ±1.0 Gyr. 

While the measurement of Vto is affected by signifi- 
cant observational errors (~ 0.05 mag), the colo r of the 
MSTO is well-defined. IChabover et all (fl996ah found 
that the usage of a point on the SBG brighter than the 
MSTO and 0.05 mag redde r (Vbto) provides more pre- 
cise relative ages than Vto- IChabover et al. ( 1996a]) pro- 
vide a conversion between Mv{BTO) and V, I data for 
a grid of five metallicities. We choose the conversion 
for [Fe/H] = — 1.5 because it is closest to the metallic- 
ity of NGC 121. In our data, we measured uibto,555 = 
22.45 ±0.02 at (m 555 -m sli ) BTO = 0.64 ±0.005 mag. To 
convert ttibto, 555 to the absolute magnitude Mbto, 555 
we use the distance modulus derived above. This yields 
Mbto,555 = 3.49 ±0.1 mag. W ith the modified cal- 
ibration by Uohnson et al.1 (|1999f ) we obtain an age of 
11.50 ± 0.5 Gyr for NGC 121. We have summarized all 
our age results in Table [3] 

4.2.3. Red Bump 

Calculating Am£gf Bump = m 555i Bump - m 5b5 HB we 
find —0.13 ± 0.05 mag. There is a general trend of 
increasing RGB bump br ightness with decreasing age 
(jAlves k Saraiedinil fl999h . assuming that a different 
"second param eter" is not affectin g the p osition of the 
RGB bump. lAlves k SaraiedinU (|1999h presented a 
^Bwmp vs - Pe/H] diagram (their Fig. 6), where the 
brightness difference between the RGB bump and the 
HB is plotted against the cluster metallicity. If we as- 
sume that the HB magnitude does not critically depend 
on age, then NGC 121 is slightly older than 10 Gyr based 
on this relation. 

These comparative results are sensitive to the abun- 



Age [Gyr] 



Method 



Reference for method 



10.9 ±0.5 
11.5 ±0.5 
10.8 ± 1.0 
11.8 ±0.5 
11.2 ±0.7 



AV HB 
M v (BTO) 
My (HB) 
Isochrones 
Isochrones 



10.5 ± 0.5 Isochrones 



Walker (19921 (ZW) 
^jhabovei^et^il^ (1996a) 
Chabover et al. (1996b) 
Pietrinferni et al. f 2004') 
Girardi, "private communication" 
Girardi et al. (2000) 
Potter et al. (2007) 



"All derived ages are listed along with the method applied. In 
all cases, we adopted the ZW84 metallicity scale. 

dance of a-elements in NGC 121 as compared to Galac- 
tic globular clusters. Many nearby Galactic globular 
cluste rs are enhanced in a- elements relative to the solar 
value ([Johnson et al.ll2004l ). As mentioned earlier, spec- 
troscopic results for NGC 121 indicate that this cluster is 
similarly enhanced in a elements, which means that the 
relative ages should not be affected as long as we con- 
fine ourselves to the comparison of globular clusters with 
similar a-element ratios. The fact that NGC 121 does 
not follow the trend of reduced [a/Fe] ratios observed in 
other nearby dwarf galaxies facilitates both our relative 
age determinations. 

4.2.4. Relative Age of NGC 121 

In Table 2] we compare the relative age of NGC 121 
with those for a sample of Galactic globular clusters. 
While this comparison sample is located in the Galactic 
halo, some objects may have formed outside the Galaxy 
and might have been subsequently captured or accreted. 
We list the clusters by their identification in column (1). 
The [Fe/H] values are given in column (2) in the scale by 



HB 

xn 



and 



iZinn k West! (fl98l . Column (3) shows t he AV n 
column (4) the ages obtained by using the Walker "(1992) 
calibration. Finally, column (5) gives the relative age dif- 
ference of these clusters to NGC 121 S(t) w For A V£? 
we adopted the values from iDe Angeli et alj (|2005l ). un- 
less differently stated (see footnotes of Table 11} . 

The clusters are listed in order of increasing AV^q 
and are divided into two groups. The first group shows 
nine "pure" GCs with similar metallicities as NGC 121. 
Among these nine clusters is N GC 2808 for which mul- 
tiple MSTOs have been found (|Piotto et aLll2007t). Its 
AVtS-P value, derived prior to the study by (|Piotto et al.1 
2007), is comparatively small. Even though those nine 
clusters have all similar metallicities, the spread in 
AV^j and therefore in age is quite large: NGC 1262 
shows the lowest AV^q = 3.24 and is similar in age 
to NGC 121, while NGC 6656 has AV^J 1 = 3.55, which 
makes it -3 Gyr older than NGC 121. 

The second group includes a subset of Galactic halo 
clusters that appear to be significantly younger than 
the average of the Galactic globu l ar cluster popula- 
tion ( e.g., iRosenberg eFail 1 19991 IVandenBergl [2000; 



iSalaris k Weiss! 12002ft . Some members of this group 
are listed in the last part of Table 2] NGC 362 and 
NGC 288 are known to be a second parameter clus- 
ter pair of differe nt ages, as reflected i n their differ- 
ent A V£? (e.g.. [Fusi Pecci et al l 119961: ICat elan et all 
[2001: iBellazzini et al.ll2001h . As NGC 362 has a sim- 
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ilar 



AV^f and [Fe/H] as NGC 121, it should there- 
fore be of a similarly young age. Other members 
of this group of y oung halo globulars are IC 4499 
(|Ferraro et al.l ll995). Ruprecht 106 (assumed to be 3- 
5 Gyr younger than the bulk of the Galactic globu- 
lars with similar metallicitics; 



: iDa Costa et a 


.1 119921: 


(Buonanno et al.l 


1995a). 


b), and Pal 14 (Saraiedini 



Terza n 7 ( Buonanno et aljl995 
[1991 . 

Among the theories that try to explain the existence 
of these young objects is the model according to which 
they are intergalactic clusters captured by the MW 
(|Buonanno et al.1 Il995b[ h or clusters formed during in- 
teractions between the MW and the M agellanic Clouds 
assum in g that they a re on bound orbits (|Fusi Pecci et al.l 
Il995h . IZinnl (|1993l ) argued that the apparent young 
halo globular clusters formed in dwarf galaxies that later 
merged with the MW. Hence the Galactic globular clus- 
ters are assumed to be a mixture of objects that formed 
with the MW itself (old halo group) and others accreted 
from destroyed dwarf satellites (young halo clusters) (see 
also Mackey & Gilmore 2004). At least six globular clus- 
ters are believed to be associated with the Sagittarius 
dwarf galaxy (e.g., ICarraro et al.|[2007l and references 
therein), providing support for the accretion scenario. 
Note, however, that Sagittarius is contributing both old 
(M54, Ter 8, Arp 2) and "young" (Ter 7, Pal 12, Whit- 
ing 1) globular clusters to the MW. Similarly, the only 
other Galactic dSph galaxy known to contain globular 
clust ers, Fornax, would contribu te both kinds of globu- 
lars (|Buonanno et al.lll9 98, 1999) if it were t o merge with 
our G alaxy. This also holds for the LMC (|01sen et al.l 
119981 see also discussion in Grebel, Gallagher, & Har- 
beck 2003). 

If we take all the ages determined in this paper into 
account, we find that NGC 121 is consistently 2-3 Gyr 
younger than the oldest Galactic g lobular clusters (ab- 
solut e age ~ 13 Gyr according to iKrauss fc Chaboverl 
|2003|) and LMC globular clusters. The age offset re- 
mains when comparing NGC 121 to old Galactic globular 
clusters in the same metallicity range (see Tab. |4j upper 
panel). We also show that NGC 121 is not as young as 
the youngest Galactic and Sagittarius globular clusters, 
some of which are ~ 2 Gyr younger than NGC 121. 

5. SUMMARY AND DISCUSSION 

We derived ages for the old SMC globular cluster 
NGC 121 based on our high dynamic range HST/ACS 
photometry that extends at least three magnitudes below 
its MSTO. In order to obtain absolute ages, we applied 
three different isochrone models. These isochrone models 
yielded ages of 11.2 ± 0.7 Gyr (Padova), 11.8 ± 0.5 Gyr 
(Teramo), and 10.5 ± 0.5 Gyr (Dartmouth). We find 
the a-enhanced Dartmouth isochrones provide the clos- 
est approximation to the MS, SGB, and RGB, whereas 
the other models cannot reproduce the slope of the up- 
per RGB. High-resolution spec troscopy indicates t hat 
NGC 121 is indeed a-enhanced (| Johnson et al.ll2004h . a 
property that it shares with many of the old outer Galac- 
tic halo globulars. Given the proximity of NGC 121 to 
the SMC on the sky and its distance, its physical associ- 
ation with the SMC seems well-established. 

Our determinations of relative ages for NGC 121 are 
consistent with the results of our absolute age determi- 



nation. Relative age estimates, when converted to an ab- 
solute age scale, are 10.9 ± 0.5 Gyr (AVfif ), 10.8 ± 1.0 
(M V {HB)) and 11.5 ±0.5 Gyr (M v{BTO) ). These num- 
bers agree well with the absolute age derivations. Our 
results confirm that NGC 121 is 2-3 Gyr younger than 
the oldest MW and LMC clusters (as also found in earlier 
WFPC2 studies). 

NGC 121 is similar in a ge to the youngest globu lar clus- 
ter in the Fornax dSph ([Buonanno et al.lll999h . and to 
several of the young Galactic halo clusters. On the other 
hand, NGC 121 is not as young as some of the Sgr dwarf 
galaxy's globular clusters or the youngest Galactic glob- 
ular clusters. 

It is intriguing that the SMC - in contrast to other 
Galactic companion dwarf galaxies with globulars - does 
not contain any old classical globular clusters. But given 
the existence of only one cluster and the question of star 
cluster survival, this could be a result of the one sur- 
vivor from the SMC's epoch of globular cluster formation 
randomly sampling an initial distribution of star cluster 
ages. On the other hand, in low-mass galaxies without 
bulges, spiral density waves, and shear it is much more 
difficult to destroy globular clusters through external ef- 
fects. That this cluster is both younger than the Galactic 
mean and enhanced in a-elements may have interesting 
implications for the early development of the SMC. 

It also is intriguing that the only globular cluster in the 
SMC is not very metal-poor. The SMC must have ex- 
perienced substantial enrichment prior to the formation 
of NGC 121. In the LMC, where two main epochs of the 
formation o f populous compact star clusters have been 
found (e.g., iBertelli et al.|[l992l ). a few globular clusters 
are found that are old enough to exhibit blue HBs. In- 
terestingly, these globular clusters, which are similarl y 
old as the oldest Galactic globulars (Ol sen et al.lll 998). 
have a similar metallicity to NGC 121 ([Johnson et al.l 
l200l (e.g., NGC 1898, NGC 2019), indicating very early 
chemical enrichment. The MW also contains old classi- 
cal globular clusters (with blue HBs) that have similarly 
high metallicities as the somewhat younger NGC 121. 
Evidently, the conditions for and the efficiency of star 
formation varied in these three galaxies at early epochs. 

After NGC 121 formed there was a hiatus in surviv- 
ing stars clusters and thus possibly in cluster formation 
activity in the SMC: The second oldest SMC cluster is 
Lindsay 1 with an age of ~ 8 Gyr (Glatt et al. 2007, 
in preparation). Since then compact populous star clus- 
ters formed fairly continuously until the present day in 
the SMC (e.g., Da Costa 2002) - in contrast to both 
the LMC and the MW. In forthcoming papers on our 
ACS photometry of SMC clusters and field populations 
we will explore the evolutionary history of the SMC in 
more detail. Clearly, clues about the early star formation 
history of the SMC will have to come from its old field 
populations. 
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Institute, which is operated by the Association of Uni- 
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TABLE 4 

Comparison of globular clusters ages (vertical 
method) 



Cluster 


[Fe/H]zws4 


AV HB 


Age [Gyr] 


6(t)w 


NGC 121 


-1.46 


3.27 


10.9 





NGC 1261 


-1.32 


3.24 


10.4 


-0.5 


NGC 5272 


-1.66 


3.24 


10.8 


0.1 


NGC 2808 


-1.36 


3.25 


10.6 


-0.3 


NGC 3201 


-1.53 


3.28 


11.0 


0.1 


NGC 5904 


-1.38 


3.34 


11.4 


0.5 


NGC 6254 


-1.55 


3.37 


11.9 


1.0 


NGC 6218 


-1.40 


3.48 


12.9 


2.0 


NGC 6752 


-1.54 


3.53 


13.7 


2.8 


NGC 6656 


-1.41 


3.55 


13.7 


2.8 


Pal 12 
IG 4499 
NGC 362 

Pal 14 
NGC 288 


-0.94 
-1.75 
-1.33 
-1.65 
-1.40 


3.17 
3.25 
3.27 
3.33 
3.45 


9.45 
11 
10.7 
11.7 
12.6 


-1.3 
0.1 
-0.2 
0.8 
1.7 



"Ages were determined using the IWalkerl lfl99l 1 cali- 
bration. The metallicity o f NGC 121 was adopted from 
Da Costa & Hatzidimitriou (1998). The other results for 
NGC 121 were derived in t his Paper. Th e data for 
NGC 6656 where taken from [R osenberg ct al. (1999); for 
IGC 4499 from IFerraro et al.1 (119951 ); for Pal 12 from 
IRutledge et al.l 119971) and for Pal 14 from IFerraro et al.l 
U995I) and IFerraro etaTl HTMEl) ; ISaraiedini (I1997TT AU 
other values were taken from Dc Angcli et al. (2005). 
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